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O ' Abstract 

We prove polynomial and exponential decay at infinity of eigen- vectors of partial 
p I I differential operators related to radiation problems for time-harmonic generalized 

' Maxwell systems in an exterior domain 17 C M^, > 1, with non-smooth inho- 

mogeneous, anisotropic coefficients converging near infinity with a rate r~'^, t > 1, 
towards the identity. As a canonical application we show that the corresponding 
eigen- values do not accumulate in M \ {0} and that by means of Eidus' limiting 
absorption principle a Fredholm alternative holds true. 
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1 Introduction 

To establish a solution theory for time-harmonic boundary value problems in exterior 
domains it is now well known that Eidus' limiting absorption principle [3] is a major tool. 
For this, one crucial step is to show that there are no point eigen-values, or at least that 
possible point eigen-values do not accumulate and that the corresponding eigen-spaces 
are finite dimensional. The absence of non vanishing eigen-vectors can be proved by a 
general pattern, which was suggested by Vogelsang [201 EI] and Eidus [3] and consists of 
the following partial results: 

step 1: eigen-solutions decay polynomially 

step 2: eigen-solutions decay exponentially 

step 3: eigen-solutions have compact support 

step 4: eigen-solutions vanish 

These results are well known, for instance, for Helmholtz' equation including perturba- 
tions. See PQ El El El ESI 113 EI] and the literature cited there. In the case of time-harmonic 
Maxwell's equations, steps 1 and 2 have been shown by Eidus |5] and step 4, the unique 
continuation property, is an old result due to Leis [TOl [HI [I2]. The pattern was just re- 
cently completed in a sufficient manner by Bauer [2], who could prove the last remaining 
step 3. All these results are known for C^-coefficients with proper decay at infinity except 
of step 1, which could have been proved even for L°^-coefficients by Picard, Week and 
Witsch [I7|. 

In the paper at hand we address to the steps 1 and 2 for a generalized time-harmonic 
Maxwell problem formulated in the language of alternating differential forms. To show 
step 1 we follow closely the arguments of Picard, Week and Witsch [T7] and step 2 will be 
proved by the methods of Eidus [3] . We note that steps 3 and 4 are still open problems 
in our general case. The only known result for step 4 is the case of scalar- valued C^- 
coefficients. 

We consider an exterior domain, i.e. a connected open set with compact complement, 
Q C M^, A > 1, as a A-dimensional Riemannian manifold with compact boundary and 
the generalized time-harmonic Maxwell equations with real frequency u ^ 



6H + iujeE = -ieF, dE + ioofiH 
L*E = 



i/iG 



in dQ, 



in Q 



(1.1) 
(1.2) 



together with the corresponding radiation condition 



{-ly^ *dTA*H + E,dTAE + H decay at infinity. 



(1.3) 
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Here E, F and H, G are differential forms of rank q (g-forms) and g + 1 {{q + l)-forms), 
respectively, and d resp. 6 = (—1)'^^ * d* is the exterior differential resp. co-differential, 
the latter acting on (g + l)-forms. By * we denote as usual Hodge's star operator and 
by A the exterior product, t : dQ "-^ Q is the natural embedding of the boundary and 
L* is the pull-back of t, i.e. the tangential trace operator. We intend to model non- 
smooth, inhomogeneous and anisotropic media by linear transformations e and /i on q- 
and {q + l)-forms, respectively. 

For the sake of a short notation we introduce the pars of q-{q + l)-forms 

u:={E,H), f:={F,G) 

and for those the formal matrix operators 



M :-- 



where i? := dr A and T := 

as 






6 


, S:= 


'0 


T 


, A:= 


'e 0" 


d 





R 





/i 



Ma 



iA-^M, 



-\Y^ *R*i and write our problem (ll.ip - fll.3p more compactly 



(Ma -uj)u = f in Q, 

L*E = indn, (1.4) 

[S + l)u decays at infinity. 

For the system (11. 4p we will show polynomial and exponential decay of eigen-forms. 
For the polynomial decay we can admit L°°-coefficients e, fi, while we need C^-coefficients 
to prove exponential decay. In both cases the coefficients must converge at infinity with 
a rate r""^, t > 1, towards homogeneous and isotropic coefficients. 

The main tool to handle irregular coefficients is a decomposition lemma, which allows 
us to prove the polynomial decay of eigen-forms by reduction to the similar result known 
for the scalar Helmholtz equation. The keys to this decomposition lemma are weighted 
Hodge-Helmholtz decompositions, i.e. decompositions into irrotational and solenoidal 
forms, in the whole space case and a well known procedure to decouple the electric and 
magnetic form by discussing a second order elliptic system. To illustrate this calculation 
let us look at (11.41) in the homogeneous case. Applying Mh + OJ yields 

{M^ + uj^)u = 0. (1.5) 

By (II. 4p E is solenoidal, i.e. SE = 0, and G irrotational, i.e. dH = 0, since 66 = 0, 
dd = 0. From A = d6 + 6 d, where the Laplacian acts on each Euclidian component, we 
get the identity M^m = {6 dE,d6H) = Au and finally (II. 5p turns to the componentwise 
Helmholtz equation 

{A + uj^)u = 0. (1.6) 

The polynomial decay of eigen-forms together with an a priori estimate for the solu- 
tions corresponding to non-real frequencies is sufficient to prove a Fredholm alternative 
for (ll.4p utilizing the limiting absorption principle invented by Eidus [3j. Moreover, we 
get at most finite dimensional eigen-spaces for possible eigen-values but these can not 
accumulate in M \ {0}. 
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2 Definitions and preliminaries 

For later purpose let us fix ro > 0, such that \flG Br^, where denotes the open 
ball of radius 6 centered at the origin. We also define the exterior of the closed ball 
Bo :— IR-^ \ Bff and the sphere So, both of radius 9. 
Using the weight function 

p := (1 + r'^y^'^, r{x) := x, 

we introduce for m e No and s e R the weighted scalar Sobolev spaces 

H™(r2) := e ll^{n) -.p'O'^^e V|a| < m). 

In Vt we have a canonical global chart, the identity, and thus, Vt becomes naturally a N- 
dimensional smooth Riemannian manifold with Cartesian coordinates {xi, . . . ,xp^}. For 
alternating differential forms of rank q E Z (g-forms) we define componentwise partial 
derivatives 9"$ = {8°' ^j) dx^ , if $ = ^jdx^ (sum convention!), where / are ordered 
multi-indices of length q. Then, for m e No and s e R we define weighted Sobolev spaces 
H^'^(f2) of g- forms as well. Equipped with their natural scalar products all these spaces 
become Hilbert spaces. For m = we also utilize the notation L^''^(0) := H°'^(0) and for 
pairs of forms we introduce product spaces like 

In the special case s = we neglect the index and we have in L2'«(0) = H°'''(0) = Ho''^(0) 
the scalar product 

Jn Jn Jn Jn 

for $ = $y dx^ , \& = ^idx^ e L^'*(Q). Here A denotes Lebesgue's measure and ( • , •)q 

the pointwise scalar product on g-forms. 

By Stokes' theorem and the product rule the exerior derivative and co-derivative are 
formally skew-adjoint to each other, i.e. 

(d$,*)L2.,+i(f^) = -($,5*)L2.,(n) V($,*) e C°°'^'^+i(Q), 
which gives rise to weak definitions of d and 6. Here we denote the vector space of all 

o 

smooth g-forms with compact support in Q by C°°'^(r2). We note that still dd = 0, 55 = 
and d 5 -|- 5 d = A hold true in the weak sense. Furthermore, for s e R we introduce some 
special weighted Sobolev spaces suited for Maxwell's equations 
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Equipped with their natural graph norms these are Hilbert spaces as well. To generalize 

o o 

the homogeneous boundary condition we introduce D^(Q) as the closure of C°°'^(Q) in the 
norm of D^(Q). Utihzing Stokes' theorem we see that in fact the homogeneous boundary 

o o 

condition l*E — is generalized in D^(Q). The spaces Df(Q), A^(Q) and even Dj{Q) are 
invariant under multiplication with bounded smooth functions (p. A subscript at the 
lower left corner indicates vanishing exterior derivative resp. co-derivative, for instance, 

o o 

oDj{n) = {$ G Dj{n) : d$ = 0}. The indices loc and vox refer as usual to local 
integrability and compact supports, respectively. If the whole space Q = is under 
consideration, we omit the dependence on the domain and write simply, for example, 
oD^ := oD^(R'^). Moreover, for weighted Sobolev spaces Vj, i e R, we define 

t<s t>s 

Now let us introduce the properties of our transformations e, /i, A: 

Definition 2.1 Let r > 0. We call a transformation e r-admissible, if 

(i) e{x) is a linear transformation on q-forms for all x e Jl; 

(ii) e possesses L°°{^)- coefficients, i.e. the matrix representation of e corresponding to 
the canonical basis (and then for every chart basis) has L°° (D,)- entries; 

(iii) e is symmetric, i.e. 

V ^ e L2'«(Q) L2,.(Q) = £*) L2,«(Q) ; 

(iv) e is uniformly positive definite, i.e. 

3c> 0V$ e L2'^(n) (£$,$)L2„(a) > cml,,(^^y, 

(v) e is asymptotically the identity, i.e. e — £oId-|-£ with £o £ ^+ o-nd e — Q{r'~'^) as 
r — >■ oo. 

Moreover, for n & Nq we call £ T-C"-admissible, if e is T-admissible and 

(vi) e G €"(-8^0) with bounded derivatives, which means that the matrix representation of 
i corresponding to the canonical basis (and then for every chart basis) has C\Bro)- 
entries and all derivatives are bounded. 

We call T the order of decay of the perturbation e. 

We remark that by a transformation x :— ax, H :— f3H we may assume with loss of 
generality £0 = A*o = 1 throughout this paper. 
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Finally, we note that the multiplication operators R, T and S are related to the 
differential operators d, S, M through the following formulas: 

CoMr) = ^'ir)X, = i rX3^, = - i 3^rX 

Here {D, X) G {(d, R), {6, T), (M, S*)} and Ca,b denotes the commutator of two operators 
A, B, i.e. CD,ip{r) = D(p{r) — (p{r)D, where (p is a smooth function on M. Furthermore, 
5" denotes the Fourier transformation on g-forms in (componentwise in Euclidean 
coordinates), which is an unitary mapping on L^''^. 

3 A decomposition lemma 

The following decomposition lemma is essential and allows us to transfer results known 
from Helmholtz' equation to Maxwell's equations without any further regularity assump- 
tions. To use results from Week and Witsch [23j we set 

I:= {n + N/2, 1-n- N/2 : n G No}. 

Lemma 3.1 Let A be r-admissible with order of decay r > 0. Furthermore, let K be a 
compact subset of C\ {0}, u E K , t, s E with < s G M \ I and t < s < t + t. Let 
6 > tq and ip := r]{r/6), where rj G C°°(R) is a cut-off function supported in [l,oo) and 
constantly equal to 1 in [2,oo). Moreover, let u G D^(f2) x A^^^(r2) satisfy 

{M^-uj)u=:f 

Then, 

f :=-iipAf + {Cm,^ - iuj^k)u G L^'^'^+i 

and by decomposing 

according to ^23, Theorem 4] 

/:= Ad+-M/sGoA^XoDfi 
holds. Then u can be decomposed into 

U= {I — Lp)u + 'Ud,5 + + Ua 

and there exist generic constants c > 0, which are independent of u, f or u, such that 
(i) (1 - ip)u G Dl^^{n) X A^+i(fi) and for all t G M 

1(1 -'^)^llD?(n)xA?+i(n) ^ c(||/||L2,.,,+i(f^) + ||m||l2^..,+1(^)); 
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(ii) u^,s := --(/d,5 + /s) e X Af 1 and 

ll^d,5|lDlxAl+i ^ c||/||l2,9.9+i; 

(iii) := 3^-i((l + r2)-i(l - ir,S)3^/) e H^'^'^+i n (oA« x qD^+i) anc/ 

(iv) UA := «5,d- e H,2'^'^+i n (oA? x oDf+^) and 

||liA||H?'«.9+i — c(||ua||l?'«'«+i + ll'^9'llHi'«'«+0 

for all t <t, where 

us,d:= --(Ad- Mipu) e Hj'^'^+i n (oA? x oDf 



These forms solve 

{M + iuj)(pu^ f, {M + iu;)us,d^ f, {M + iuj)uA ^ {I - iuj)u3r, 
(M + 1K = / 

and 

(A + uj^)ua = -(1 + oj'^)us + (1 - ic^)/. 
Moreover, the estimates 

||/|||_2.g., + l < C||/|||_2,,,g + l, 

||/|||_2,9., + 1 < c(||/||L2.<7.g+l(^) + \\u\\^2,g.^g + l^^^), 



and 

\\{ A + uj'^)u a\\^2,<,,,+i < c(||/||L2,«,<;+i(n) + \\u\\i_2,,,,+i(^^-^) 
as well as 

||(M - i A5)m||l2,,,,+1(^) < c(|||/||L2,,,.+i(f^) + ||m||l2^,,,+1(^) + ||(M - i A,5)ma||l2,.,.+i) 
hold for all i <t and uniformly in X & K, u and f. 
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Proof Obviously, (Ma — u)u = / G L^'''^+^(f2) is equivalent to 

{M + iuA)u= - iAf e L2'5'9+1(^]). 
Since (pu G x Af^^ we have 

Mipu = ipMu + Cm,<pU = — i ouipAu — i ipAf + Cm,lpU 

and thus 

(M + iu;Vn = /GL^''?''^+\ (3.1) 

since Cm,lp = t]' {r j 6)6~^ S is compactly supported and t + t > s. We rewrite (13. ip in the 
form 

lUipU = -Mipu + fs,d+ fd,5 + fs 

and note 

us,d= --{fsA-M^u) e (D^noAf) X (Af^ noDf 1) c H,i'^'^+\ 
Wd,5 = --(/d,5 + /s)eD^x Afi 

with {pu = usA+ Md,5 and by regularity, e.g. [9^ Lemma 4.2(i)]. (For s < N/2 we even 
have /s = 0.) Moreover, w^^d solves 

Mus,d = M(pu - Mud,5 = -iujus,d+ f5,d+ -Mfs, 

CO 

i.e. (M + iuj)u5,d = f e oAj x oD1+^. Now, to define (M + l)-^ by the Fourier 
transformation we put 

Then, e L^'^-^+i as well as Ju^^ G l^''^'^^^ are implied by 3^/ G L^'^-^+i. Hence, 

From / G L^'^'''"'"^ we get by definition 3^f G H'*'^''^'*"^. The components of J'uy arise from 
those of 3^f by multiplication with bounded C°°-functions. Thus, also 

G H^'«'''+^ 

follows; see e.g. Wloka [25| p. 71, Lemma 3.2]. Again, by definition, ujr g L^'^'''"'"^ and we 
obtain the estimate 

||Mg-|||_2,9,<j+l < c||/ II |_2, 5,9+1. 

Since M3^~^ = iS^^^rS we compute 

(1 + r^)3'{M + l)uj = (1 + ir^)(l - ir5) J/ = (1 + r^S^)3'f. 
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SF — and dG — imply TUF — and R3^G — 0, respectively, and therefore, using 
RT + TR^l 



TR 
RT 



holds. Hence, 

J{M + l)uy = i.e. (M + l)uy = /. 
Besides, we have G (D« fl qA^) x (A^+i n oD^^) and thus, 

again by regularity. Considering 

UA = «5,d- e H,^'«'^+i n (oA? X oDfi) 

we calculate 

(M + ia;)'UA = (1 — iuj)uj. 
Once more by regularity we even obtain ua G Ht^'«''^+^ n (oA? X oDf+^) and we compute 

(A + uj^)ua = (M -iuj)(M + iuj)uA = (1 - ia;)(M - ic^)iigr 
= -(l+a;>:r+(l-ia;)/. 

Finally, we achieve the asserted estimates from the regularity result and the continuity of 
the projections in L^'^ onto oDf, qA^ resp. mentioning that 



\\Mfs\\ < c\\fs\\ 

holds in any norm since S^''^"'"^ is finite dimensional and M linear. 



□ 



4 Polynomial decay 

First, we need a trivial but useful technical lemma. 

Lemma 4.1 For all t, t G M with i < t and all i!} > there exist constants c,9 > 0, such 
that 

holds for all ijj e 

Proof For sufficient large 6' > we get from t — t < 

Thus, lim (1 + ^^)*~* = completes the proof. □ 
Our decomposition lemma implies: 
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Theorem 4.2 Let A be t -admissible with r > 1. Moreover, let I C M± be a closed 
interval and u I as well as 1/2 < s G M \ I. // 

^2 > 2 

is a solution of Maxwell's equation 

(MA-a;)«=:/eL^^+i(fi), 

then u e Dg_]^(r2) x AglJ(r2) and there exist constants c,9 > independent of u, f oroo, 
such that 

Proof Let t > -1/2 and u E D?(f2) x A^+\f2) with t < s-1. Without loss of generahty 
we may assume t + l<s<t + r. Otherwise, we replace t and s by := t + ka and 
Sk := t + 1 + {k + l)a < s, k = 0, . . . , with a := (r — l)/2 > and obtain the assertions 
after finitely many a-steps. 

Decomposing u by Lemma 13.11 we get solutions ua E H^''^''^^^ of Helmholtz' equation 



{A + uj')uA =: /a e oA^ X oDf ^ 

A componentwise application of [211 Lemma 5] yields ua E H^^f^^ and with a constant 
c > independent of ua, /a or u we have 

||Wa||h2''J,9+1 < c(||/a|||_2,9,9 + 1 + ||Ma||l2''J^<?+i) . 



Furthermore, from Lemma [3. II we have u E A^_^(i7) x Al_^(il) and the estimate 

< c(||/a|Il2.9,.+i + \\f\\i_fi^i+\n) + 

< c(||/|||_2,,,,+i(^) + ||M||L2^9.g+i(f^)), 

where we assumed without loss of generality r < 2. Since s — r < s — 1 the assertion 
follows now by Lemma 14. 1[ □ 



Remark 4.3 Let the assumptions of Theorem \4.S\ be satisfied. If f E L^''^''^+-'^(f2) for all 
s eM., we get 

u E Dun) X Ai^\n) 

for all s E M.. This holds, for instance, if f is exponentially decaying or even compactly 
supported. 
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5 Exponential decay 



Using the 'partial integration' technique introduced by Eidus [5J for the classical Maxwell 
equations we will prove: 

Theorem 5.1 Let cj G M \ {0} and A be t-C"^- admissible with r > 1. // 

> 2 > 2 

is a solution of Maxwell's equation 
for all t G M, then 

e'^^ue (D'^in) X A'^+l(^])) n h2'«'''+1(5,„+i) 

holds for all t G M. The assertion holds in particular if f is compactly supported. 

Proof The idea of the proof is to estimate the exponential series. For this, we need some 
technical preliminaries. For all s G M we have 

(M + iujA)u = - iA/ G L2'9'«+1(1]) n H2'«'9+1(5,J. 

Hence, by Theorem i^] and Remark SSI belongs to Df(fi) x A^+^(J1) for all s G M. 
Therefore, inner regularity, e.g. a combination of a cutting technique together with [21 
Lemma 4.2(i)], yields u G H^'«'''+^(5,.o+e) for all s G M and all positive 6. Thus, 



(M + iw)M = - i A/ - icuAu =: f G L^'9'''+i(f]) n HI''''''+\R 



ro+e) 



for all ^ > and all s G M. Consequently, applying 6, d and (M — iu) to the latter 
equation we can compute in Bro 



and 



iu6E = 6F, iujdH = dG 



(5.1) 



By defining □ := A — we note 
Therefore, we achieve 



'6d 


" 




'dS 


" 





d5 


, □ = 





5d 



+ □ = A. 



(A + u^)u = (M - iw - -□)/ =: /, 
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which is the equation we intend to work with. Now, we multiply this equation and all 
forms by r™' with some m G R and indicate the resulting forms by an index m. We note 
that all occurring forms are well defined elements of H^''''^"'"^(i?ro+6') for all ^ > and all 
m G M. Using Lemma [B. II we obtain in 5^ 



r, mim + 2 — N) s m„ 

{A + U^ + '-)Um - 2— dr Um 

fm = {M — iu □)/m + (— CM.r™ H Cci,r"^) f 

1 1 Tfl Tfl ~\~ 1 ~ 

-i (iM + a; + -□ (i + sw + sw swS^))f„ 



(5.2) 



With rj from Lemma ISTT] we define the cut-off function ifg := ri{r — 6 + 1) for all 6' > ro + 1. 
Then, 

supp Lfe C Be, supp Vy^e, supp(l - i^e) n Bg C Bg n Bg+i. 

Without loss of generality let any form be real-valued. We multiply f l5.2p by Lpgr^Um with 
p G M resp. ipgrdrUm and integrate over M^. We achieve 



r 



)Um - 2^ dr Um, feT^'Um) 



L2,9,9 + l 



+ m{m + 2 - N)ipgrP \Um\q^q+l - rrKpgr^ dr \Um\q^q+l dX 



(5.3) 



resp. 



Or Urn) 



m{m + 2- N). 



m 



)Ur 



- 2— dr Um, ^er dr Um) 



L2,9,9 + l 



AUm^fgr dr Um + -^^^r dr \Um\q,q+l 

jYi -|- 2 — TV 

+ V^e dr \Um\l,q+l " 2mV9e | drUm\l^q+i d\. 



(5.4) 



By partial integration we get from (15. 3p resp. (15.41) 



<fgrP[{uj + ^)| 



u 



|«|=i 



Ur 



) dX 



< I {L, <^er^M^)L2..,9+i I + cm{0 + 1)^^+^ 



(5.5) 



Polynomial and Exponential Decay of Eigen-Forms 



13 



resp. 



^ II-, 

l«l=l 



(5.6) 



+ 4m| drUmllg^i) d\ 

< I {fm, ^Per dr Mm)L2.g,9+l | + Cm^{9 + 1)2™+!^ 

where 

7m,p,jv := m{m + p) + p{p + N - 2)/2, 7^,^ := m(m + 2 - A^)(A^ - 2) 

and c is a generic constant independent of m and 6. Now, we multiply (15.51) for p = by 
N — 2 and add (15.51) and (15. 6p in a suitable way. We obtain 

r m(N — 2)2 

< c( I (/m, V30Um)L2,<7,g+i I + I (/m, dr U^) L2.^:<i+^ | + 771^(6* + 1)2™+1 



which yields 

/ ^e\uJl,+idX<{N -2f I ^Ji'^\u,^\l^^^dX + cm\e + lf^+^ 
JRN ' J^N r r 



To take care of the right hand side we prove two lemmas. 



(5.7) 



Lemma 5.2 For all p G M there exists a constant c > 0, such that for all ip G H'^^^i^ro+i) 
and all 9 > + 1 

o 

Proof By a cutting-technique we may assume ^ ^l/2i^ro+i)- Hence, it suffices to 

o 

prove the assertion only for ip G C°°'''{Bro+i) by continuity. But for those the result is 
simply shown by several partial integrations. □ 



Lemma 5.3 Let p,p,p G M with p > p and p > 2. There are a constant c > and a 
non-negative function x tending to zero at infinity, such that for a// m G M and a > 0, 
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6 > Tq + 1 the following estimates hold: 
(i) |(/^,<^,r^'OL2„„+i| <x(e)(m(e + l)2-+P+ / ^erP\Ul^^,d\ 



|a|<2 



(ii) I (/^, y.,rPM^)L2.„+i I < x(^) (m(^ + 1)2™+P + V / ^eA 9" /^Ij f/A 

■^1^"' |a|<l ^ 

(iii) |(/^,<^era,n^)L2„„+i| <x(0)(m2(e + l)2-+i+ V / (^gr'^l 9" /^Ij dA 

|«|<1 |c«l<2 

Proof By several partial integrations we remove all derivatives from kum resp. jm and 
use the decay of A resp. the integrability of /„. This yields (i). In (ii) there are no 
longer second derivatives allowed on u^- Hence, we have to insert a a into the estimate 
since we do not demand any decay properties of the derivatives of A. To prove (iii) we 
have to handle one challenging term, i.e. 

The part of the first summand causing the biggest difficulties is 

because we can not integrate by parts anymore since E is only twice weakly different iable. 
But with flS.ip and Lemma IB. II (ii) we can substitute 

Tfl 1 Tfl TD 

6E^ = -TE„. - -{5 - -T)eF^ -{6- -T)iE^. 
r Lo r r 

Thus, the most challenging term reads now 

{6eE,^,iPprdr6iE^)i2,g-i = \ I '^pr dr\5eEm\l_i dX 
and can be handled easily by two partial integrations. □ 
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We proceed with the proof of the theorem by further estimating fl5.7p using Lemma 
Oim and (iii). 



|a|<l l«l=2' 



cm'(e + l)2'"+^ 



Utihzing Lemma [5 .21 we estimate the second derivatives of Um by the first ones and Au^ 
Then we substitute Am^ with fl5.2p and get for sufficient large 6 



V / i^e\d''um\l,q+i dX 



4 



l«l<2 

Now we insert this estimate into fl5.8D and obtain 



/ ^9\Um\lg+l dX 
a|<l 

+ J2 [ ^eAl + ^)\d'^fm\l,+,dX)+cm\e + lf 



\a\<2 

For all p G M f l5.5p yields the estimate 



2 

m 



(5.9) 



|a|=l ■ 

< c / ^erP(l + dX + I (/„, (/?er%r«)L2.?..+i | + cm{e + 1)2"^+?, 
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such that by Lemma 15.31 (ii) for sufficient small a and large 6 

I 1 ■''K^ 



|a|=l ' 



+ cV / ^er^\d''Ul^^,d\ + cmie + l) 



2m+p 



\a\<l ' 

follows. Now, we plug the latter estimate for p = and p = — 4 into ( 15.91) obtaining 



<X(^)(/ ^,(1 + —)\Um\lg+, d\ 

+ V/ ^em^rP\d''Ulg_,,dX)+cm\9+l) 



2m+l 

where we assume without loss of generality p < p. Therefore, for sufficient large 6 

'fe\Um\l^q+l d\ 

+ V / ¥.emV|9"/„|2^+if/A)+cm5(^+l)2™+\ 

/ dX 

J Be 

<>ci9){ I "^Wmllq^.dX+Y^ I m'r^\d^f^\l^^,dX)+cm\e + l) 

J Be \a\<2'^ ^0 

Setting k := 2m we finally get for all A; G M and sufficient large 9 

r''Hl,q+l dX 



I.e. 



2m+l 



Be 







|a|<2 ' 


Be 



(5.10) 
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We are ready to prove the exponential decay. Let t G R+ and 6 large enough, such 
that f lS.lOp is holds. Then, we have for all natural numbers 8 < i^i < K2 



t f 

k=Ki ^ 



k=K, \a\<2 o^sW k=K, <(5(e+i))fe(e+i) 



^2 ,fc /. 



|a|<2 



Now, let 9 be so large, such that x{9)t^ < 1/2. Then, 

^2 ,k p Ki-l ^,^g 



6t{9+l) 



l^. J Ba 1 ri J Bo 



k=Ki " A:=0 



Since the right hand side is independent of K2 we obtain by the monotone convergence 
theorem for K2 — t- 00 

tr \ \2 1 \ 

e "A < 00, 



Bo 



i.e. e^^ u G L^''''^+^(i7) for all t G M. The differential equation yields 

e*" Mu G L2'«'«+l(^]), e*''(feE, d/iif) G L2'''-i'''+i(fi,„+e) 
for all t G M and all 6' > 0. Consequently, by 
Me*"M = e*'' MM + te*"^M, 

5e e*" E = e*" 6eE + t e*" TeE, d fi e''' H = e''' d fiH + t e*" 

and inner regularity (see the beginning of the proof) we achieve e^^ u G H^''^''^^^{Bro+e) for 
alH G M and all 6' > 0. Repeating this argument yields finally the same assertion for H^. □ 
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6 Solution theory for time-harmonic Maxwell equa- 
tions 

As a canonical application we intend to present a solution theory for the radiation problem 
( II ■4p . We will follow in close lines the first part of [17J. It can be seen easily that 

are skew-adjoint to each other. Consequently, for 0-admissible coefficients A 

M : D{M) := □"(fi) X A«+^(fi) C aI^''''^\^) ^ aI^''''^\^), Mu := M^u 
is self-adjoint. Here, := L'^''^''^~^^{Q) equipped with the scalar product 

{U,v) (AM,f)L2.<j,9+i(Q). 

Furthermore, we will denote the kernel and the image of M by A^(M) and /(M), respec- 
tively. 

Definition 6.1 Let a; G C \ R and f G l.f^'^^''^'^'^ [Q) . Then u solves Max(A, cu, /), if and 
only if u E i5(M) and {M\ — u)u = f . 

The self-adjointness of M yields the unique solvability of Max(A, w, /) for non-real 
frequencies w G C \ M and right hand sides / G L^''''^"''^(f2). We denote the continuous 
solution operator by C^. Since the spectrum of M is contained in the real axis we expect 
from well known facts about Helmholtz' equation that we have to work in weighted L^- 
spaces and utilize radiating solutions to get a proper solution theory for real frequencies. 

Definition 6.2 Let u G M\{0} and f G lf;^l''^^\Q) . Then u solves Max(A, /), if and 
only if 

M G DVi (fi) X A''+^ (S + 1)m e L^'^^+^(fi), (Ma - uj)u = f. 

< 2 < 2 > 2 

The second constraint will be called 'Maxwell incoming radiation condition' or simply 
'radiation condition'. 

We will establish a solution theory using Eidus' limiting absorption principle. The key 
tool for the application of this principle is an a priori estimate, which ensures the uniform 
continuity of C^j operating in proper Hilbert spaces even up to the real axis. 
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6.1 An a priori estimate 

Lemma 6.3 Let / C M \ {0} be a compact interval and s, —t > 1/2. Furthermore, let 
A be T-admissible with r > 1. Then, there exist constants c, ^ > and at > —1/2, 
such that for all u G C+, which means that u has got non negative imaginary part, with 
u;^ = + iaA, X e I, a e {0, 1] and f e as well as u := f the estimate 

holds true. 



Proof Without loss of generality let s G (1/2, 1). We note (1/2, 1) fll = 0. Decomposing 
/ and u using Lemma [3.11 with s = s, t = we have /s = since s < N/2 and obtain 
Ma € H^'''*+^ satisfying 

(A + u^)uA = -(1 + u;^)uj + (1 - iu)f =: /a G L2'«'^+\ 

The self-adjointness of A : H^'^-^+i c L^'^.^+i ^ i2,Q,g+i yields (A + w^)-!/^ = ua- 
Applying \24l Lemma 7] , which is a well known a priori estimate for the scalar Helmholtz 
equation in M^; see also Ikebe and Saito ^ or Vogelsang [20l section 2], componentwise 
to Ua and using Lemma \3A] with M(e~^^^MA) = e~'^^(M — iXS)uA we get the estimate 

||Ma||l2>9.9 + 1 + \\iM — iAS')'UA II I 2,5,9+1 
^ s— 1 



< C^||Ma||l2.9.9+1 + ^ II 5"(e Ua)|||_2,9,9 + 1 

H=i 

< c||/a||l2>9.9+1 < C(||/|||_2,,,5+l(f^) + ||m||l2^,,,+ 1(^^). 



(6.1) 



which holds uniformly in ua, /a and A. But actually we are interested in estimating the 
term ||(M — ic<;S')MA||| 2,, ,9+1. This needs an additional argument. The standard resolvent 
estimate yields 

t'"|A| ||'Ua||l2>9.9+1 < II /a II L2. 9, 9 + 1 (6.2) 

and we have |u; + A| > |A| since | Rea;| > |A|-\/2/2 and lu G C+. Thus, by (16. 2p and 

, u^-X^ iaX 
oj - X = — = - 

uj + X uj + X 

we achieve uniformly in u 



\{M — ia;S')MA||, 2,9,9+1 < ||(M — i AS')ua||, 2,9,9+1 + c\uj — A|||ua||| 2,9,9+1 

c 

< (M - iAS')MA ,2,9,9+1 + tttII/a||l2.9,9+i. 
'"-s-i A 
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A combination of tlie latter estimate with fl6.ip and Lemma 13.11 yield 

i^iD?(n)xA?+i{n) + II (^^ - 1^*^)^112^7+1 (H) 

< c^||ma||l2,9,9+i + II (M — ia;S')MA 11^2,9,9+1 + II/IIl2'«'''+1(q) + II'"IIl2'9'«+1(q)) 

< C(||/||L2,9,9 + I(f^) + ||M||L2^9,9+l(f^)) 

uniformly in u, f and u. Using (M — iujS)u = — iuj{S + 1 + A.)u + / we finally arrive at 

ll«llD?{C)xA|+i(n) + l|('^+ l)"llLtV+'{n) - c(||/IIl2*^+1(Q) + Mtl-^-^^+^n))- 

Because of the monotone dependence of the weighted L^-norms on the weights we may 
assume t near to —1/2 and s near to 1/2, such that 1 < s — t < r holds. Then Lemma 
14.11 completes the proof. □ 



6.2 Fredholm theory using the Umiting absorption principle 

We prove three more technical lemmas. 

Lemma 6.4 Let a, /3 G M with < a < (3 and \ Q G B^- Moreover, for some t G 
let E G D?(fi), H G At^\n) and ^ G C\[a,(3],C)- Then with 

: [0,/3] C 

(p{s) ds 



max{a, cr} 



and ^ := <f o r , \[f :=■?/' o r 

{(^RE , H)^2,q+if^B^r\B!i) = dE, if)L2,9+i(nns^i) + (^-E, 5-ff)L2,9(nnB^)- 

Proof Assume E G C°°''^{n), H G C'^''^+\n). With 7 := / (p{s) ds we have V'Ip.a] = 7, 
= and ip G C^((«,/3),C) with = — By Stokes' theorem we compute 

(*dE,/7)L2,9+i(cni?^) + (^'E,5i/)L2,9(nnB^) 
= 7(dE,if)L2,9+i(nnBc) +7(^,'^^)L2,9(onSc) 

= 7 / C(E A + (<l>i?E, //)L2,9+i(^„nB,) 
-7/ LUEA*H)+ij{f3) f l;{EA*H), 

J Sa J 
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where Lg : Sg ^ M denotes the natural embedding. With the help of mollifiers we get 

1 ° ° 

the desired formula for all H G A^^ Since C°°''^(f2) is dense in D^(fi) the assertion 

holds as stated. □ 

By the same approximation technique we obtain the rule of partial integration for 
weighted forms. 

Lemma 6.5 Let t, s e R and E e D?(f2), H G A9+^(f2) as well as (fg := l-r]{- /O), 
$61 := ifg or with some 6 > tq. Then 

($,dE,i7)L2,,+i(f,) + {<l>gE,SH)t^,in) = -{v>'eir)RE,H)^,^^^,^^^ 

holds. Additionally, ift + s>0 then 

{dE,H)L2.,+^n) + {E,6H)l2,,^^) = 0. 

Remark 6.6 If u e Dl{n) x Al^\n) and v G D«(fi) x A^+^fi) with t + s>0 then 

Since M is self-adjoint we have the well known Hodge-Helmholtz decomposition. 
Lemma 6.7 

|_2,g,g+l(^) = aL2'5'9+1(^]) = iV(M) ©a7(M) 

= (oD^(^^) X oA'^\^)) ©A A"^(5A9+i(fi) X dD'?(^])) 

Here ©a denotes the orthogonal sum in \L'^''''''^^{Q) and the closures are taken in the 
respective L"^ -spaces. 

Another essential ingredient of the solution theory generating convergence in the lim- 
iting absorption argument is the so called Maxwell local compactness property MLCP, i.e. 
the embeddings 

have to be compact for alH < s and all q. For a detailed analysis of this property of dQ 
we refer to [22] , [Uj and [131 El US US] as well as the papers cited there. 

Definition 6.8 For w G C \ {0} we define 

P := {cj G C \ {0} : Max(A, 0) has a nontrivial solution.}, 
:= {u : u is a solution of Max(A, 0).}. 

We remark P C M \ {0} or A^"^ = A^(M - a;) = {0} if a; G C \ M. We are ready to 
prove the main result of this section. 



22 



Dirk Pauly 



Theorem 6.9 Let w G M \ {0} and A be r-admissihle with r > 1. 

(i) Eig en- solutions decay polynomially, i.e. for all t 

Additionally, let Q possess the MLCP. Then 

(ii) is finite dimensional; 

(iii) P has no accumulation point m R \ {0}; 

(iv) for every f G L^t^^^(^) there exists a solution u of the problem Max(A, a;, /), if 
and only if fJ-\N^, i.e. 

WveN^ (A/,t;)L2.,.,+i(n) =0. (6.3) 
The solution can be chosen, such that U-LaN^^j, i.e. 

WveN^ (AM,t;)L2,,.,+i(n) = 0. (6.4) 

By this condition u is uniquely determined and the solution operator C^^ f := u is 

continuous. More precisely, maps L2''?'9+1(^]) niVj-* to (D^(fi) x A^+^(^])) HiVj-'^ 
continuously for all s, —t > 1/2. Here we denote the orthogonality corresponding to 
the \i.'^''^''^~^^{Q) -scalar product by ±a. 

Proof We follow the proof of [T71 Theorem 2.10]. To show (i), i.e. the polynomial decay 
of any eigen-solution u, we only have to prove 

u G L^t'i''\n) 

> 2 

because of Theorem 14. 2[ Remark 14. 3[ the equation Mu = — i coAu and the inclusions 



d Di{n) c oD'^+^n), 6Ai+\n) c oA'^(^^)- (6.5) 
Using the second part of the radiation condition we obtain some t > —1/2, such that 



lim \\RE + ifll, 2,9+1, o „p s < oo 



holds true. We calculate 



+ + 2Re(<l>i?E,/7)L2,9+i(B^^nB^) 
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with $ := p^*. Lemma [6.41 the differential equation and the symmetry of e, /i yield 

{(^RE,H)^2,q+i(^Br^f^Bp) = (^di?,if)L2,?+i(nnB^) + (^-E, 5/f)L2.<j(ons^) 

= -iu (^'/iif,if)L2,g+i(QnB^) + iw {'^E,eE)^_2,,^^^B^) . 

eK m 

Thus, Re{<^RE,H) 

1^,1+1 (B^^^nBis) = and hence, by means of the monotone convergence 
theorem H G i.f''^~^^{Q) follows for /3 — )■ oo. Finally we get E e r'^i(n) using the first 

^ 2 

part of the radiation condition. 

If (ii) or (iii) would be wrong then there would exist a sequence of eigen-values 
{ui)i^f^ C M \ {0} tending to u and a sequence of eigen-forms (■u^)^gN C N^^,, such that 
{u£)iQ^ is an ortho-normal system with respect to the AL^'^'^'''^(f^)-scalar product. As an 
ortho-normal system ('u^)^gN converges in L^''^''^^^(f2) weakly to zero. Moreover, by the 
differential equation {ui)em is bounded in 

Hence, from the MLCP we can extract a subsequence {uT^e)e^fi, where tt : N — t- N is strictly 
monotone, converging for all t < in L^''^''^^^(f2) to 0. The latter is due to the weak 
convergence. For 1 < s G M\I Theorem 14.21 vields uniformly in (u,r£)^eN and (u;^^)^^^ the 
estimate 



1 — (Au7r£, U^^)L2,<7.g+i(n) < c||M^^||L2,q,q+l(f^) < c||M7r£||^2,Y+i 

^11 II 2 ^— >oo ^ 



which is a contradiction. 

We prove (iv): First of all (16. 3p is necessary since we get for all eigen-forms v G 
by their polynomial decay and Remark 16.61 

(A/,t;)L2.9.g+i(n) = i((M + iwA)M,f)^2,,,,+i(f^) = - i (m, (M + i ujA)v ) |_2,,.,+i = 0. 

=0 

To show existence we now use Eidus' principle of limiting absorption. For that purpose 
let / G l^'l'^~^\n) with §3^. Moreover, let {ce)£e'N be a positive sequence tending to zero 

and {fe)eeN C. L^''^''^"'"-^(f2) with some s > 1/2 be a sequence satisfying fi^^N^^, such that 
{fe)ef^fn converges to / in L^'''''"^^(i7) as i tends to infinity. Defining non-real frequencies in 
the upper half plane ue G C+ \ R with uj = u"^ + i aeu and — )■ w we obtain L^-solutions 

:= C^, fe G D{M) = D«(l]) x A^+\n) 
solving the Maxwell problem Max(A, u;^, fe), i.e. 

(Ma - ue)ue = U (6.6) 
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Applying the orthogonal Hodge-Helmholtz decomposition we decompose ui and ac- 
cording to Lemma 16.71 



ue = < + ul, f, = /f + // e N{M) ©A /(M). 

Therefore, 

(Ma - iOeWe - fe = + /f e Ar(M) n /(M) = {0}. 

As orthogonal projections the sequences of forms {f^)een, {fi)em converge in L'^''^''^~^^{Q). 
Hence, so does {uf)^^^. Let us assume the boundedness of (mD^sn or 

Vt < -l/2 3c> 0V£ G N ||m^||l2.<7..+1(j^) < c (6.7) 

for a moment. At the end of the proof we will show by contradiction that in fact (16.71) 
holds. Let t' be such a t with (16.71) . Then, (-uD^gN is bounded in 1.'^;'^''^'^^ (Q) and by the 
differential equation, i.e. 

(Ma - cue)ui = //, 

and by ([HSD even in (D^,(^]) ^£~loA^,(^])) x ( A^+^(fi) ^/i-^D^+^(^])) . Hence, the MLCP 
yields for an arbitrary i < t' a. subsequence ('U^^)^^^ converging in L?''^'''^^(i7) and even 

o ^ 

in D?(i7) X A?^ (i7) by the differential equation. Therefore, the entire sequence (mtt 

o 

converges in D?(r2) x A?^ (i7) to, let us say, 

u e D?(fi) X A?+^(fi), 

which solves 

(Ma - cu)u = f. 



With the polynomial decay of eigen-solutions and Remark 16.61 we compute for all eigen- 
forms V G A^,^ and all £ G N 

= (A/rf,t;)L2,g,9+i(f7) = -i (m^^, (M + iu;^£A)t;)^_2,,,,+i(j^) = (cJ„(? - w) (Am^£, t;)L2,,,<,+i(n). 

Consequently, (Am^^, t>)L2, 9,9+1 (n) = 0. Since ( ■ , At>)L2,9,9+i(Q) is a continuous linear func- 
tional on L?'''^^^(i7) for all v G we obtain 

VwgA^^ (AM,t;)L2,9.9+i(n) = 0. (6.8) 



Now, we pick some t < —1/2. Then, we get by Lemma 16.31 constants t > — 1/2 and 
c,6 > 0, such that by the monotone convergence theorem 



\\u 



D?(Q)xA?+\n) + < c(||/||L2.9.9+i(f^) + \\u\\L2,,,,+i{nnBg)) 



holds. Therefore, u G D'' _i(fi) x A'^~^^i{Q) and u satisfies the radiation condition, i.e. 



{S + l)u G In other words, u solves Max(A, w, /). 
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We note that this proves the principle of hmiting absorption to hold. The choice 
fi := / for all £ G N yields the existence of a solution of Max(A, /) and this one is 
unique because of fl6.8p . 

Moreover, for —t, s > 1/2 the solution operator C^^ maps Ds{CS) to It{Cuj), where 

D^ic^) := L^.''.''+l(^]) n h{c^) := (D?(^]) X Ar\n)) n iv^, 

continuously. This follows by the closed graph theorem because Da{C^) and It{C^) are 
Hilbert spaces by the polynomial decay of eigen-solutions and C^j is closed. The latter 
assertion is a consequence of Lemma 16.31 and the monotone convergence theorem. 

Finally, it remains to contradict the contrary assumption to (16. 7p . To this end, let 

o 

t < -1/2 and (M^)£eN C D^(^]) x (^) be a sequence with ||M£|||_2,5,g+i^^j — )■ oo. Defining 
the normalized forms 

ue := \\ut\\ll,^,+,^^^ut, fe := \\ui\\-l,,,+,^^^fi 
we have ||'U£|||_2,q,9+i^j^^ = 1 for all £ G N and ||/^ 1^2,9,9+1(5^) —J- 0. Moreover, the equation 

(Ma - ue)ue = fe 

holds. Following the arguments above we obtain a subsequence {un£)£m converging in 
l_2,<?,g+ij^^-j with i < t towards u G N;^^, which solves Max(A, w, 0). Hence, u = and 

Lemma [63] yields constants c,9 > independent of a-^e, fne or u.„e, such that 

1 = i'^TT^IlL^'i-'J+Vn) — ^{ ll/'r^llL^''''«+\n) + ll'^7r£||L2,9,9+i(f7nBe) ) 
' "•^ ' v ' 

holds true; a contradiction. □ 

The polynomial decay of eigen-solutions proved above and Theorem 15.11 yield: 

Corollary 6.10 Let w G M \ {0} and A be t-C'^ -admissible with r > 1. Then, any 
eig en- solution u G A''^^ decays exponentially, i.e. 

holds for all t eM.. 

Remark 6.11 The polynomial resp. exponential decay of eigen-solutions holds for arbi- 
trary exterior domains Q, i.e. Q does not need to have the MLCP. 

Remark 6.12 // the medium is homogeneous and isotropic in the exterior of some ball, 
i.e. supp A U (M^ C Be for some 6 > 0, then 

u = in Be 
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for all uj E M.\ {0} and u G A^^^. Because in this case u solves componentwise Helmholtz' 
equation (A + ijj'^)u = in Bq and therefore, by Rellich's estimate 118] or 112, p. 59], 
must vanish in Bg. If the principle of unique continuation holds for our Maxwell system 
under consideration then even 

= {0}. 

Using the a priori estimate of the hmiting absorption principle and some indirect 
arguments followed by the 'trivial' decomposition of L^'''(n) from [13 Lemma 4.1] we 
are able to prove stronger estimates for the solution operator £^ as the ones given in 
Theorem 16.91 (iv). We only note the results here. For this, let Q possess the MLCP and 
let s, —t > 1/2 as well as be a compact subset of C+ \ {0}. Furthermore, let A be 
r- admissible with some r > 1. 

Lemma 6.13 There exist constants c,9 > and some t > —1/2, such that 

W^^ f\\Dl(n)xA''+\n) + f\\q,i.i+^(^^) < c(||/||L2.9,g+i(j^) + || fWl^^i^i+^innBe)) 

holds for allooeK and f e L2'9''2+i(fi) n N^^. 

Corollary 6.14 Let K Pi P = ^. Then there exist constants c > and t > —1/2, such 
that the estimate 

i /iD?(n)xA?+^(Q) + II ("^ + ^) fWlfi'i+^n) - '^i/iL^«'''+^(n) 

holds true for all u & K and f G L^''''^~'"^(f2). In particular, the solution operator C^^ 

o 

mapping L^''''^+^(f2) to Dj(fi) x Aj"*" is equi- continuous with respect to to E K. 
Theorem 6.15 Let K D P = ^. Then, the mapping 

w I — > 

is uniformly continuous. Here we denote the set of hounded linear operators from some 
normed space X to some normed space Y by B{X,Y). 



A Appendix: classical equations 

We want to point out briefly, which classical equations are covered by our generalized 
approach. Since the relation between the differential form calculus and the classical 
vector calculus is very well known we directly translate our equations (11. 4p . i.e. in the 
longer version fll.ip - fll.3l) . into terms of vector analysis. 



Polynomial and Exponential Decay of Eigen-Forms 



27 



g = 0: E, F are scalar functions and H, G are vector fields. We have the equations of 
linear acoustics with Dirichlet boundary condition in first order form 

divH + i ueE = — i eF, gradE" + i ujiH = — i jiG in il, 

E = ondfl 

and ^-H + E, E^ + H decay at infinity, where ^{x) := x/r{x). Using the differential 
equation for H, i.e. 

-iuj{C - H + E) ^^■gTadE-iujE+ 

=drE 

we see that we get Sommerfeld's classical radiation condition, i.e. {dr — iuj)E decays 
at infinity. 

q — N — 1: E, F are vector fields and H, G are scalar functions. We have the equations 
of linear acoustics with Neumann boundary condition in first order form 

gradi? + i ueE — — i eF, divE + i lo/jlH — — i //G in Jl, 

u-E^O on on 

and + E, ^ ■ E + H decay at infinity. Here u is the unit normal vector pointing 
outwards. Now, using the differential equation for E we get again Sommerfeld's ra- 
diation condition. Note that by the differential equation the homogeneous boundary 
condition is equivalent to the inhomogeneous Neumann boundary condition 

di,H — —\ujv ■ e{E — f) —\v ■ F on 

q — N: E, F are scalar functions and H, G vanish. We have trivial equations 

uE — —F in Q, E decays at infinity. 

q — 1, N = E, F, H, G arc all vector fields. We have the classical Maxwell equations 
with homogeneous electric boundary condition, i.e. dfl is a perfect conductor, 

— curlif + i ueE = — i eF, curli5 + iujiiH = — i fiG in f2, 

V X E = Q OTO. do. 



and — ^ X H+E, ^ x E+H decay at infinity. The latter are the classical Silver-Miiller 
radiation conditions for Maxwell's equations. 
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B Appendix: some calculations 

Let tr and sw denote the formal trace and the formal swapping operator, respectively, i.e. 



tr 



A a 
b B 



A + B, sw 





'A 


a 




'B 


a 


sw 


h 


B 




h 


A 



We note 



TR 
RT 



RT 
TR 



and formally 1 - -5^ = sw -5^ and tr S"^ ^tisw S'^ ^ RT + TR^ I. Moreover, 

SM 



'Td 


" 




'SR 


" 





R5 


, MS = 





dT 



and formally tr 5'M = Td+RS, tr MS — dT + 6R. By a straight forward computation 
we obtain: 



Lemma B.l Let m e M. 



(i) 




^ —1 m 


(ii) 




= — Rr^ , Cg^r"^ — — ^ 


(iii) 




= '^{5R + Td 


(iv) 




^"^{MS + SM '"^ + 
r r 


(v) 




772 

= — (swM5 + sw5M- 
r 


(vi) 




777 

= (trM,5 + tr5M 
r 



m 

-( 

r 



m + 2-N, 

r ' 



m + 1 



r 

m + 1 



sw 5* )r 



Looking at the two formulas for C/^^^r"' we get the nice relation 

dT + 6R + Td+R6 = tr{SM + MS) = {N - l)/r + 2dr . 
Here dr is meant componentwise in Cartesian coordinates. 
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